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Abstract 

This paper is concerned with stochastic processes that model multiple (or it- 
erated) scattering in classical mechanical systems of billiard type, defined below. 
From a given (deterministic) system of billiard type, a random process with tran- 
sition probabilities operator P is introduced by assuming that some of the dy- 
namical variables are random with prescribed probability distributions. Of par- 
ticular interest are systems with weak scattering, which are associated to para- 
<— j i metric families of operators Pj lt depending on a geometric or mechanical param- 

Q-T eter h, that approaches the identity as h goes to 0. It is shown that (P/, - I)/h 

converges for small h to a second order elliptic differential operator L on com- 
pactly supported functions and that the Markov chain process associated to P/, 
converges to a diffusion with infinitesimal generator £. Both P/, and L are self- 
adjoint (densely) defined on the space L 2 (H,77) of square-integrable functions 
over the (lower) half- space H in U m , where rj is a stationary measure. This mea- 
sure's density is either (post- collision) Maxwell-Boltzmann distribution or Knud- 
sen cosine law, and the random processes with infinitesimal generator L respec- 
tively correspond to what we call MB diffusion and (generalized) Legendre diffu- 
sion. Concrete examples of simple mechanical systems are given and illustrated 
by numerically simulating the random processes. 

1 Introduction 

The purpose of this section is to explain informally the nature of the results that will be stated 
in detail and greater generality in the course of the paper. 

A type of idealized multi- scattering experiment is depicted in Figure ITTTI The figure repre- 
sents the flight of a molecule between two parallel solid plates. At each collision, the molecule 
impinges on the surface of a plate with a velocity v and, after interacting with the surface in 
some way (which will be explicitly described by a mechanical model), it scatters away with 
a post-collision velocity V. The single scattering event v <-> V, for some specified molecule- 
surface interaction model, is given by a random map in the following sense. Let D-0 denote the 
half-space of vectors v = {y\, V2, v^) with negative third component. It is convenient to also 
regard the scattered velocity V as a vector in D-0 by identifying vectors that differ only by the 
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sign of their third component. A scattering event is then represented by a map from H into the 
space of probability measures on H, which we call for now the scattering map; the probability 
measure associated to v is the law of the random variable V. Thus the scattering map encodes 
the "microscopic" mechanism of molecule-surface interaction in the form of a random map, 
whose iteration provides the information about velocities needed to determine the sample 
trajectories of the molecule. 

The mechanical- geometric interaction models specifying the scattering map will be lim- 
ited in this paper to what we call a mechanical system of billiard type. Essentially it is a con- 
servative classical mechanical system without "soft" potentials. Interactions between moving 
masses (comprising the "wall sub-system" and the "molecule sub-system," using the language 
of |7]) are billiard-like elastic collisions. 




Figure 1.1: An idealized molecular flight between two solid plates, as an example of a multi-scattering 
experiment. We refer to v and V, respectively, as the pre- and post-scattering velocities at a 
collision event, and regard V as a random function of v, as explained in the text. 



An example of a very simple interaction mechanism (in dimension 2) is shown in Figure 
13.11 That figure can be thought to represent a choice of wall "microstructure." In addition 
to a choice of mechanical system representing the wall microstructure, the specification of 
a scattering map requires fixing a statistical kinetic state of this microstructure prior to each 
collision. For the example of Figure 13.11 one possible specification may be as follows: (1) 
the precise position on the horizontal axis (the dashed line of the figure) where the molecule 
enters the zone of interaction is random, uniformly distributed over the period of the peri- 
odic surface contour; (2) at the same time that the molecule crosses the dashed line (which 
arbitrarily sets the boundary of the interaction zone), the position and velocity of the up- 
and-down moving wall are chosen randomly from prescribed probability distributions. The 
most natural are the uniform distribution (over a small interval) for the position, and a one- 
dimensional normal distribution for the velocity, with mean zero and constant variance. (The 
variance specifies the wall temperature, as will be seen.) In fact, one general assumption of 
the main theorems essentially amounts to the constituent masses of the wall sub -system hav- 
ing velocities which are normally distributed and in a state of equilibrium (specifically en- 
ergy equipartition is assumed). In this respect, a random-mechanical model of "heat bath" 
is explicitly given. Once the random pre-collision conditions are set, the mechanical system 
describing the interaction evolves deterministically to produce V. Note that a single collision 
event may consist of several "billiard collisions" at the "microscopic level." 

Having specified a scattering map (by the choices of a mechanical system and the constant 
pre-collision statistical state of the wall), a random dynamical system on H is defined, which 
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can then be studied from the perspective of the theory of Markov chains on general state 
spaces (fT3l). 

Clearly, one can equally well envision a multiple scattering set-up similar to the one de- 
picted in Figure [O] but inside a cylindrical channel or a spherical container rather than two 
parallel plates; or, more generally, inside a solid container of irregular shape, in which case a 
"random change of frames" operator must be composed with the scattering operator to ac- 
count for the changing orientation of the inner surface of the container at different collision 
points. (See j8] . This is not needed in the case of plates, cylinders, and spheres.) We like to 
think of this general set-up as defining a random billiard system, an idea that is nicely illus- 
trated by |5][6], f° r example. 

We are particularly interested in situations that exhibit weak scattering, in the sense that 
the probability distribution of V is concentrated near v or, what amounts to the same thing, 
the scattering is nearly specular. Our systems will typically depend on a parameter la that 
indicates the strength of the scattering, and we are mainly concerned with the limit of the 
velocity (Markov) process as h approaches 0. This will lead to novel types of diffusion pro- 
cesses canonically associated with the underlying mechanical systems. We call h the flatness 
parameter for reasons that will soon become obvious. 

For the systems of billiard type considered here (introduced in Section|2), the essential in- 
formation concerning their mechanical and probabilistic definition is contained in two linear 
maps: C and A on U m+k , where k is the number of "hidden" independent variables (whose 
statistical states are prescribed by the model) and m is the number of "observed" variables 
(say, the 3 velocity coordinates of the molecule in the situation of Figure ll.l) . This m is also 
the general dimension of H. Vectors in D-0 will be written v = (v\,...,v m ). The maps C and 
A are non-negative definite and Hermitian; C is a covariance matrix for the hidden velocities 
and, by the equipartition assumption, it is a scalar multiple of an orthogonal projection, while 
A contains (in the limit h -* 0) information about the system geometry and mass distribution. 

The first observation (which is studied in much greater generality in [7]) is that the resulting 
Markov chain on velocity space H has canonical stationary distributions given by what we 
refer to as the post-collision Maxwell-Boltzmann distribution of velocities. (The term "post- 
collision" is used to distinguish it from the more commonly known distribution of velocities 
sampled at random times, not necessarily on the wall surface.) This velocity distribution has 
the form 



where a 2 = Tr(CA)/Tr(A A ), A A is the restriction of A to the subspace of "hidden velocities," 
and dV denotes Euclidean volume element. The scattering map can be represented as a (very 
generally) self-adjoint operator oni 2 (H,/i) ofnorm 1, which we indicate by where fa is the 
flatness parameter. We denote the density of \i by g := dp/dV. The term \v m \ in g equals the 
speed times the cosine of the angle between v and the normal to the scattering surface; this 
cosine factor is often referred to in the applied literature as the Knudsen cosine law. ((2]) 

Let (H) denote the space of compactly supported smooth functions on the half-space. 
A first order differential operator can be defined on this space using C and A as follows: 



(1.1) 




(1.2) 



(DO) ( v) := y/l [A 1 / 2 ( v m grad y ®-<S> m { v) v) + Tr (CA) 1/2 <£ OT ( v) e] 
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where e is the coordinate vector (0, . . . , 0, 1 ) and the subindex m in <t> m indicates partial deriva- 
tive with respect to v m . We now define a second order differential operator on (H) by 

L:=-D*D, 

where D* indicates the adjoint of D with respect to the natural inner product on the pre- 
Hilbert space of smooth, compactly supported square integrable vector fields with the Maxwell- 
Boltzmann measure /i. We refer to L as the MB-Laplacian of the mechanical-probabilistic 
model, D as the MB-gradient, and -D* the MB-divergence. 

The central result of the paper is that, as h approaches (and under commonly satisfied 
further conditions to be spelled out later), the Markov chain process with transition probabil- 
ities operator P; ( converges to a diffusion process on H whose infinitesimal generator is the 
MB-Laplacian £ . The resulting process, which we call MB-diffusion, is illustrated in a number 
of concrete examples in the paper. 

-v 2 
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Figure 1.2: Sample trajectory of the MB-diffusion process whose infinitesimal generator is given by 
(£<£)( p) = -2vi<&i + [-4i/2 + (l + v \ ) / v 2] $2 + L'f't'll -2l>i V2<$i2 + (1+ v\)^22 obtained 
by simple Euler approximation. We have used a time interval of length 50, initial condition 
(0, -1) and number of steps 50000. The parameters chosen here are not related to those of 
Figure nOl so the axis scales are not comparable. 

The MB-diffusion can be expressed as an Ito stochastic differential equation 

dV t = Z( V t ) dt+b(V t )dB t , 
where B t is m-dimensional Brownian motion (restricted to H), Z( v) is the vector field 

Z( v) := -2Av + (qJq) [(Ay, v) +Tr(CA)] e 
and b( v) is the linear map 

b(v)u:= v m A l l 2 u - ( A 1/2 v, u) e m + Tr (CA) 1/2 u m e. 
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A sample path of an MB -diffusion in dimension 2 is shown in Figure [L2l 

It is interesting to note that, in dimension 1, the operator L reduces to the (up to a constant) 
Laguerre differential operator (in this case on functions defined on the interval (-00, 0)): 

1 Id/ d®\ 

where g( v) = a~ 2 i/exp (- i/ 2 /2er 2 ) is the Maxwell-Boltzmann density and A is the scalar equal 
to the (in this case 1-by-l) matrix denoted above by A. 




Figure 1.3: Sample path for the Legendre diffusion in dimension 2 with eigenvalues (of A) X\ = 2.5 and 
A2 = 1. The starting point is (0,0), the time length is 5, and the number of steps is 50000. 
Note that diffusion is faster along the horizontal axis as Ai > A2 ■ The stationary distribution 
for this process is the normalized Lebesgue measure; therefore, long trajectories fill the disc 
evenly. 

It has been implicitly assumed above that the number of independent "hidden" velocity 
components k is positive. The case k = is somewhat different and has special interest. Now, 
particle speed does not change after scattering (we refer to this case as random elastic scat- 
tering) and both the Markov chain process and the diffusion approximation can be restricted 
to a unit hemisphere in H. Alternatively, by orthogonal projection from the hemisphere to the 
unit open ball in dimension m - 1, we can consider these processes taking place on the ball. 
Now £ is completely determined by A and has the form 

(£0)(i/) = 2f ^((l-M 2 )*,), 

where the A, are the eigenvalues of A and <1> is a compactly supported smooth function on the 
unit ball. (We have chosen coordinates adapted to the eigenvectors of A.) 
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The associated diffusion process, written as an Ito stochastic differential equation, has the 
form 

dV t = -4AV t dt+ [2(l - \V t \ 2 ) A] 1/2 dB t , 

where B t is (m - 1) -dimensional Brownian motion restricted to the unit ball. The operator L 
in this case naturally generalizes the standard Legendre differential operator on the unit inter- 
val (-1,1); we call the stochastic process on the higher dimensional balls Legendre diffusions. 
The stationary measure turns out to be standard Lebesgue measure on the ball, so Legendre 
diffusions have the interesting property that sample paths fill the ball uniformly with proba- 
bility 1. A sample path of a Legendre diffusion process is illustrated in Figure fOl One can also 
think of the Legendre diffusion as a special case of the MB-diffusion in the sense explained in 
Proposition[5]and in the remarks immediately after this proposition. 

The relationship between the scattering operators P and the above differential operators 
of Sturm-Liouville type suggests that one should be able fruitfully to investigate the spectral 
theory of P based on an analysis of L and a spectral perturbation approach. A very simple 
observation in this regard is indicated in [9], while [10 1 discusses the spectral gap of P (which 
can often be shown to be a Hilbert-Schmidt operator) in very special cases. We hope to turn 
to a more detailed analysis of the spectrum of P in a future study. 

2 Mechanical systems of billiard type 

We introduce in this section the main definitions and basic facts concerning classical me- 
chanical systems of billiard type and their derived random systems. Component masses of 
a given mechanical model interact via elastic scattering that admit a billiard representation. 
Particular attention is given to weak scattering, in which reflection in this billiard representa- 
tion is nearly specular. A sequence of random scattering events comprises a Markov chain on 
velocity space whose transition probabilities operator, in the case of weak scattering, is close 
to the identity. 

2.1 Deterministic scattering events 

The reader may like to keep in view the examples of Figures[2T|and[3T|while reading the be- 
low definitions. For our purposes, a system of billiard type is a mechanical system defined 
by the geodesic motion of a point particle in a Riemannian manifold of dimension n with 
piecewise smooth boundary. Upon hitting the boundary, trajectories reflect back into the in- 
terior of the manifold according to ordinary specular reflection and continue along a geodesic 
path. Except for passing references to more general situations, the configuration manifolds 
of the systems considered in this paper are Euclidean. More specifically, we consider ( n + 1 ) - 
dimensional submanifolds of T" x [R with boundary, for some n, with a Riemannian metric 
which will have constant coefficients with respect to the standard coordinate system. The 
boundary, assumed to be the graph of a piecewise smooth function, and the metric coeffi- 
cients are the distinguishing features of each model. The periodicity implied by the torus fac- 
tor is a more restrictive condition than really needed, but these manifolds are a natural first 
step and describe a variety of situations of special interest. Thus the configuration manifold 
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M is assumed to have the form 

M = { (x, x n+ i ) e T" x R : x„+i > F(x)} 

where F is a piecewise smooth function. 

The Riemannian metric on M is specified by the kinetic energy quadratic form, which de- 
pends on the distribution of masses in the system. The Euclidean condition means, in ef- 
fect, that the kinetic energy form becomes, after a linear coordinate change, the standard 
dot-product norm restricted to (the tangent bundle of) M, while the torus component of M 
has the form T" = n" =1 (U/ <z,-Z) for positive constants a\,...,a n . 

A scattering event is defined as follows. Let c be an arbitrary constant satisfying F(x) < c for 
all x e J n . The submanifold x n+ \ = c will be called the reference plane. We identify the tangent 
space to M at any point on the reference plane with U n+1 and denote by H" +1 the lower-half 
space in R" +1 , which consists of tangent vectors whose (n+ l)st coordinate is negative. 

Definition 1 (Deterministic scattering event). A scattering event is an iteration of the corre- 
spondence (x,v) (x', V), where x, x' lie on the reference plane and (x', V) is the end state of 
a billiard trajectory that begins atx with velocity v and ends atx' with velocity V . By reflecting 
V on the reference plane, we may when convenient regard both v and V as vectors in H" +1 . 

Notice that the map describing a scattering event is indeed well defined, at least for almost 
all (x, v), by Poincare recurrence. If Igrad^Fl is uniformly small over x e T", a condition that 
is assumed in the main theorems, trajectories cannot get trapped. 

The iteration of the scattering event map introduced in DefinitionQ] as well as its associ- 
ated random maps described in Subsection l2.2l acquires greater significance in the context of 
random billiard systems as in |7], but the various concrete examples given later in this paper 
(the simplest of which appears in Subsection l2.3> should provide enough motivation. 

2.2 Random systems and weak scattering 

The random scattering set-up defined here is a special case of the one considered in (7). 
Briefly, the main idea is that some of the variables involved in a deterministic scattering event, 
as defined above, are taken to be random. The scattering map then becomes a random func- 
tion of the initial state of the system. The resulting random system can model a variety of 
physical situations; we refer to |7] for more details on the physical interpretation. 

The notation 7(X) will be used below to designate the space of probability measures on a 
measurable space X. We start with a deterministic scattering system with configuration man- 
ifold McT"xR and boundary function F : J" -> R. Recall that M is defined by the inequality 
x n+ \ > F(x), x e T". The deterministic scattering map is then (x, v) >-* (x', V), where x,x' 
lie on the reference plane x n+ \ = c (recall that c is an essentially arbitrary value that specifies 
the reference plane); v and V lie in the lower-half space H" +1 , and V is the reflection on the 
reference plane of the velocity of the billiard trajectory with initial state (x, v) at the moment 
the trajectory returns to the reference plane. A scattering event can consist of several billiard 
collisions. 

Choose c' such that sup^. |F(x)| <c'<c and define 

M c , := {(x,x„+i) € M: x n+ i > c'} = T" x (c',oo). 
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Let k < n be a non-negative integer and write M c i = J k x j n ~ k x (c',oo). Accordingly, decom- 
pose the tangent space to M at any point on the reference plane as H" +1 = R k x H™ _ . Fix 
a probability measure p on [R* and set m = n - k+ 1. By a random initial state with observ- 
able component v e we mean a state of the form (x, c, w, v), where x e T" is a uniformly 
distributed random variable, c is the value defining the reference plane, and w is a random 
variable taking values in R k with probability measure p. To this random initial state we can as- 
sociate a probability measure v v € T(1H!! 7 ) as follows: Consider the trajectory of the system of 
billiard type having random initial state (x, c, w, v), and let V be the component in H" 7 of the 
final velocity of this trajectory, reflected back into U™, at the moment the trajectory returns 
to the reference plane. Then V is a random variable and v v is by definition its probability 
measure. We refer to v v as the return probability distribution associated to the random initial 
state having observable component v. 

Definition 2 (Random scattering event). Let\i be a probability measure onU k and give J" the 
uniform probability measure, denoted A. Then the random scattering event associated to the 
system of billiard type and these fixed measures is defined by the map 



where v v is the return probability associated to the random initial state with observable com- 
ponent v. 

The probability measure p on U k typically will be assumed to have zero mean, non-singular 
covariance matrix of finite norm, and finite moments of order 3, when not assumed more con- 
cretely to be Gaussian. The uniform distribution on T" is, by definition, the unique translation 
invariant probability measure. 

Definition 3 (Scattering operator P). Let Co(H™) denote the space of compactly supported 
continuous functions on the lower half-space. For any given $ e Cq(H™), define 



We call P the scattering operator of the system for a random initial state specified by p and the 
uniform distribution on the torus. 

Operators similar to our P naturally arise in kinetic theory of gases and are used to spec- 
ify boundary conditions for the Boltzmann equation. See, e.g., Q][TT]. Typically, the models 
of gas surface interaction used in the Boltzmann equation literature are phenomenological, 
such as the Maxwell model (| 1 1, Equation 1.10.20), and are not derived from explicit mechan- 
ical interaction models as we are interested in doing here. 

From the definitions it follows that P and v„ are related by 



where the expression in the middle denotes the expectation of the random variable <$(V) 
given the initial condition v. 



(P®)(v):= f f ®(V(x,c,v,w))dA(x)du(w) 
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Based on the examples given throughout the paper, we can expect P and v h> v v generally 
to have good measurability properties, due to the deterministic map from which the random 
process is defined being typically piecewise smooth. In our general theorems it will be im- 
plicitly assumed that v v v (A) is Borel measurable for all Borel measurable subsets A of 
H™. Billiard maps are typically not continuous; see (3) for basic facts on billiard dynamics (in 
dimension 2) . 

The following additional assumption turns out to be convenient and not too restrictive. 

Definition 4 (Symmetric M, F). The configuration manifold M of a random scattering process 
or, equivalently, the function F defining it, will be called symmetric ifF(o+ u) = F(o- u) for 
all ueU" and some choice of origin o in T" . 

2.3 Example: collision of a rigid body and flat floor 

A simple example will help to clarify and motivate some of the above definitions. (More rep- 
resentative examples will be introduced later.) Consider the 2-dimensional system of Figure 
12.11 It consists of a rigid body in dimension 2 of constant density and mass m that moves in 
the half-plane set by a hard straight floor. There are no potentials (e.g., gravity). The body and 
floor surfaces are assumed to be physically smooth, in the sense that there is no change in the 
component of the linear momentum tangential to the floor after a collision. The motion of 
the center of mass can then be restricted to the dashed line of Figure l2.1l due to conservation 
of the horizontal component of the linear momentum. 







X 






F(9) 





Figure 2.1: Collision of a rigid body and a flat floor in dimension 2. The shape of the body is encoded in 
the function F(8). For a disc, F is constant. Some of the main results of the paper apply to 

shapes for which h := supg (F '(0)) is small. 

Figure F2\2l shows the description of the same example explicitly as a system of billiard type. 
Let B represent the body at a fixed position, with its center of mass at the origin. Define the 
second moment of the position vector by I 2 := Area -1 f B \b\ 2 dA(b), where A is the area mea- 
sure. Set coordinates X\ = 6 and xz = x/ 1, where 8 is the angle of rotation and x is the height 
of the center of mass of the body at a given configuration in U 2 . Then the configuration man- 
ifold of the system is the region M= {(xi,Xz) e Tx R : F(xi) < IX2}, equipped with the kinetic 
energy metric K=k (xf + xf ), where k is a positive constant. We model the collision between 
the body and the floor by a linear map C : T X M -> T X M, where x is a boundary (collision) point 
of M. Under the assumption of energy conservation and time reversibility, C is an orthogonal 
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involution; the assumption of physically smooth contact is interpreted as Cu = u for every 
nonzero vector u tangent to the boundary at x. As C cannot be the identity map, it must be 
standard Euclidean reflection, whence the system is of billiard type. 

The 3-dimensional version of this example is similarly described, the function F now being 
defined on the special orthogonal group SO(3). The kinetic energy metric on M c SO(3) x U 
is no longer Euclidean and naturally involves the body's moment of inertia. 



s 




2tc 

Figure 2.2: The billiard representation associated to the mechanical system of Figure [2TTI Here M c 
T x IR. (Dashed lines indicate periodic conditions.) The reference plane is indicated by the 
top horizontal line and the wavy ground is the graph of F. We define the pre-scattering angle 
f and initial position s e [0, 2n\ The outgoing angle is E. 

One way in which this deterministic system can be turned into an example of a random 
system is by regarding the initial angle 8, at the moment the center of mass of the body crosses 
a reference plane, to be random with the uniform distribution over the interval [0,2;r]. In 
other words, suppose that the exact orientation of the body in space at a given moment prior 
to collision is completely unknown. As the magnitude of the velocity of the billiard particle 
(that is, of the moving point particle of Figure 12.21 is invariant throughout the process due 
to energy conservation, we may consider the return probability as being supported on the 
half-circle in Hp , which we identify with the interval of angles [0, n] . 

Thus the probability distribution of the return angle S given <f is the measure: 

l r 2n 

U~v ( (U) = — lu(~(s,0)ds 
In Jo 

where U is a measurable subset of the interval [0,n] (a set of scattered angles) and ly is the 
indicator function of U. Similarly, given a continuous function <D on [0,n], 

1 r 2j[ r It 

(PO)(0 = -/ <S>(Z(s,0)ds= *(S)v f (H). 
2n Jo Jo 

The probability distributions of the velocity of the center of mass and the angular velocity (ex- 
pressed in terms of xz and ±\, respectively) are obtained by taking the push-forward of un- 
der the maps u |<f|cos(w) and u ^ |£|sin(«), respectively. Note that approaches weakly 
the delta measure 8^ supported on £ when the body becomes more and more round (hence 
the reflecting line in Figure l2T2l becomes more and more straight). In this case, P approaches 
the identity operator. 
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2.4 Further notations 



All the examples of deterministic systems given in this paper can be turned into random scat- 
tering systems in various ways. The most natural choices of random variables fall within the 
scope of the following discussion, in which the definition of a random scattering event is re- 
stated in a more convenient form. Let the tangent space to M at any point (x, c) of the refer- 
ence plane decompose in the following two different ways: 

H" +1 = R k x H m = R" x (-oo,0). 

where m = n - k+ 1. Accordingly, any given £ € T X M has components relative to these two 
decompositions defined by 

Let {e\, . . . , e n+ \) be the standard basis of IR" +1 and e:= e n+ \ the last basis vector. So 

where the inner product represented by the angle brackets is the standard dot product. The 
component £ Y of the final velocity of the billiard trajectory is the quantity of interest produced 
by the scattering event. The component £ A of the initial velocity is assumed random with a 
probability distribution /i. 




Figure 2.3: The random collision process. The initial position r e J n x {c} on the reference plane is 
chosen randomly with the uniform distribution. The outgoing velocity in TM is f and 
the vertical projection V := f Y is the outgoing observed velocity. We thus obtain a random 
map v i-» V. It will turn out to be convenient to express both r and V in terms of ( and the 
independent variable x e T", as indicated in the figure. 



We can now express the random scattering map v <-> V by the following algorithm, which is 
illustrated in Figure l2~3l 
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Definition 5 (Random scattering algorithm). In the notation introduced above, the random 
scattering map is defined by the following steps: 

i. Start with v e U™; 

ii. Choose a random weU k with the probability distribution p, and form <f := w + v; 

Hi. Choose a random r = (7, c) , where 7 e T" is uniformly distributed, and let ( r, be the 
initial state for the billiard trajectory; 

iv. Leti,' e H" +1 be the velocity of the billiard trajectory at the time of its return to r n+ \ = c; 

v. SetV:=(£') y eH™. 

For billiard surfaces (i.e., the graph of F) that are relatively flat, the typical collision process 
comprises a single collision. It makes sense in this case to introduce the independent variable 
x as indicated in Figure 1231 and allow both 7 and to be functions of x and f . For xel", let 
n{x) be the unit normal vector to the graph of F at (x,F(x)). Note that 

e-grad x F 

nix) = - - . 

Vl+lgrad^P 

As a geometric measure of the strength of scattering we introduce the following parameter. 

Definition 6 (Flatness parameter). The quantity h := sup^j,, Igrad^Fp will be referred to as 
the flatness parameter of the system defined by F. 



3 Stationary measures and general properties of P 

The basic properties of P are described in this section. These are mostly special cases of re- 
sults from 1 7 1 , which we add here for easy reference. Proofs are much simpler in our present 
setting and are sketched here. 



3.1 Stationary measures 

It will be assumed in much of the rest of the paper that the probability distribution fi for the 
velocity component w in U k , when k > 0, is Gaussian: 

e -\\w\ 2 la 2 

dp( w) = —j-dV( w) 

(2na 2 ) 1 

occasionally referring to a 2 as the temperature (of the "hidden state"). Let A be the transla- 
tion invariant probability measure on T". The standard volume element in open subsets of 
Euclidean space will be written dV, or dV k if we wish to be explicit about the dimension. 

Recall that the deterministic scattering map T is defined as the return billiard flow map on 
the phase space restricted to the reference plane, that is, T" x H" +1 = (T" x IR fc ) x H" ! , where 
we are factoring out the observable velocity component W from the hidden component at 
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temperature a 2 , which is given the probability distribution v := A® ^ in the sense described 
in the previous section. Denote by n the natural projection n : T" x U x HP -» ¥\™ 

When k = (no velocity components among the hidden variables), the scattering interac- 
tion does not change the magnitude of the velocity in HP; thus one may restrict the state to 
the unit hemisphere S" !_1 in HP. Let dco(u) represent the Euclidean volume element (mea- 
sure) on the hemisphere at the unit vector u. When necessary we indicate the dimension of 
the unit hemisphere as da> m ~ l . Observe that 

dV m (v) = c\v\ m - x d(jj m -\v/\v\) d\v\, 

where c is m times the ratio of the volume of the unit m-ball by the volume of the unit ( m - 1 ) - 
sphere. 

The Markov operator P naturally acts on probability measures on HP as follows. With the 
notation i](f) := f f dr\, the action of P on 77 is the measure rjP such that (f]P)(f) = i](Pf), 
for every compactly supported continuous /. A probability measure r\ on U™ is said to be 
stationary for a Markov operator with state space H" 1 if rjP = 77. 

The action of P on probability measures has the following convenient expression. Given 
any probability measure 77 on HP, we can form the probability measure v®rj on (T" xR fc ) xH™, 
then act on this measure by the push-forward operation under the return map, and finally 
project the resulting probability measure back to HP. (We recall that T*(, for a given measure 
C, can be defined by its evaluation on continuous functions as (T*()(f) := T).) The 
result is r\P. 

Lemma 1 . The operation r\ rjP for rj e CP (HP ) can be expressed as 

T]P={no r)*(v®77), 

where v = A ® \i is the fixed probability on the hidden variables space J n x U k , T is the return 
map to the phase space restricted to reference plane, identified with J" x [Rr x HP, and n is the 
projection from this phase space to HP. 

Proof. The straightforward proof amounts to interpreting the definition of P given earlier in 
terms of the push-forward notation. See (7) for more details. □ 

Proposition 1. When k = 0, the measure dr){ v) = ( v, e) da>( v) defined on the unit hemisphere 
in H" ! is stationary under P. Identifying the unit hemisphere with D" !_1 = {x e RP" 1 : \x\ < 1} 
under the linear projection (x,x m ) h* x, the stationary probability is, in this case, the normal- 
ized Lebesgue measure on Df _1 . Fork>0, the measure 

dn{v) = \{v,e)\e-\\ v \ 2 l a2 dV m {v) 

on HP is stationary under P. 

Proof. For a much more general result see [7] . We briefly show here the second claim. First 
note that the measure d( (£) := (<f,e) dX(x) dV(^) on T" x H" +1 is T-invariant. (The term 
(£, e) contains the cosine factor that appears in the canonical invariant measure of billiard sys- 
tems in general dimension.) The measure d((£) = exp(-||^| 2 /cr 2 ) drjo(^) is also T-invariant, 
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since any function of |£| is invariant under the return map T. Now, consider the decomposi- 
tion £ = (x, w, v) e T" x K k x H™, under which £ splits as £ = v® 77, where 

rf7 7 (^) = |(i;,e)|e-2l , 'l 2 /' 72 dy m (y). 

Here we have used: (£, e) = ( v, e) and the splitting of the exponential involving |£| 2 = | w\ 2 + 1 v\ 2 
as a product of exponentials in w and v. Thus 7r*£ = 77, where 7r* indicates the push-forward 
operation on probability measures. We now apply Lemma[T] noting that T*£ = £, to obtain 

tjP=(jio r)*.(v®77) =7i*T*( = 7i*( = ri, 

which is the claim. □ 

There is a significant literature in both pure mathematics and physics /engineering con- 
cerning random billiards, in which ordinary specular billiard reflection is replaced with a ran- 
dom reflection. The typical assumption is that the post-collision velocity distribution corre- 
sponds to the above Maxwellian distribution or, more simply, to the Knudsen cosine law with 
constant speed. See, for example, |2][5]|6]- 

It is natural to regard P as an operator on the Hilbert space L 2 (H'",ri) of square integrable 
functions on the observable factor with the stationary measure given in Proposition Q] The 
next (easy) proposition is proved in (7). 

Proposition 2. Suppose that the billiard system is symmetric, as defined in Definition^ and let 
77 be one of the stationary measures described in Proposition^ Then P : L 2 (H™, 77) -* L 2 (H!", rf) 
is a self-adjoint operator of norm 1. 



3.2 Example: collision between particle and moving surface 

The example discussed here is the simplest that exhibits most of the features of the general 
case. Its components are a point mass m\ and a wall that is allowed to move up and down. See 
Figure l3Tl The wall surface, which could be of any dimension n > 0, has a periodic, piecewise 
smooth contour and the up and down motion is restricted to an interval [0, ao/2]. In the 
interior of this interval the wall moves freely, bouncing off elastically at the heights and ao/2. 
Collisions between the wall and mass m\ are also elastic. 




Figure 3.1: Collision between the surface of a moving rigid wall and a point mass m\. The entire wall is 
given a certain mass mo; models with more localized mass definition are also possible. 
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Let x denote the coordinate along the (horizontal) base of the wall and let Zq be the height 
at which the base stands at any given moment relative to its lowest position. The range of 
Zq is assumed to be [0, «o/2]. The contour of the wall top surface, when Zq = 0, is described 
by a periodic function f(x) of period a\. Thus, when the base is at height Zq, that contour 
is the graph of x i-> f(x) + z®. It is convenient to allow zq to vary over the symmetric interval 
[-ao/2, ao/2] and set the wall surface function of x, zq as G(zq,x) = f(x) + \zq\, which can 
then be extended periodically over K 2 . The graph of G so extended is, up to a rescaling of 
the coordinates to be described shortly, the surface shown in Figure [3721 Periodicity of G is 
expressed by G(x+ ma\,za + nao) = G(x,zq) for integers m, n. Equivalently, we think of G as 
a function on the 2-torus. 




Figure 3.2: Billiard representation for the moving wall example. 



The coordinates of the point mass m\ are represented by (z\,zz) respectively along and 
perpendicular to the base line of the wall. Thus the state of the system at any moment is 
specified by ( zq , z\ , z 2 , vq , v\ , v 2 ) , where vq is the velocity of mo and v = ( v\ , V2 ) is the velocity 
of m\. 

An appropriate choice of coordinates makes the kinetic energy metric explicitly Euclidean. 
Set x = (xo,xi,xz), where xq := \J7noJ~m\ zq/ a\ , x\ := z\\a\, X2 ■- Z2I a\. The above function G 
in this new system becomes 

F(xq,xi) = a[ l f(aixi) + \/m]mo\xo\. 

Defining t := ^y^< then F(xq + rm,x\ + n) = F(xq,x\) for integers m, n. The configuration 
manifold of the particle-movable wall system can now be written in terms of F as 

M = {x € U 3 : x 2 > F(x , xi )}. 

The kinetic energy of the system then becomes K(x,x) = Xo||^|| 2 > where Kq = mia 2 /2. Under 
the assumption that the wall surface is physically smooth, we obtain again a system of billiard 
type in M, as depicted in Figure l372l 

A random billiard scattering process based on the above set up can now be defined as fol- 
lows. The observable state space is the set H_ of approaching velocities, consisting of the 
vectors u\ e\ + 112 e 2 , u 2 < 0. The part of the phase space of the deterministic process on which 
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the return map T is defined is T 2 x {0} x R x H 2 . Observe that T 2 = (K/tZ) x (IR/Z) has co- 
ordinate functions (xo,Xi), and the reference plane, with equation %2 = 0, is identified with 
T 2 . At the initial moment of the scattering event it is assumed that the height of the wall (xo) 
and the position of mi along a period interval of the wall contour (xi) are random uniformly 
distributed over the respective ranges. Thus the initial position on T 2 is a random variable 
distributed according to the normalized translation-invariant measure. 

Also at the initial moment of the scattering event the velocity of the wall is assumed to be a 
Gaussian random variable with zero mean and variance a\. That is, the initial derivative w := 

xo of xo is normally distributed with mean and variance a 2 = j^-^z- Thus the probability 

distribution for w is given by the measure [i such that dfi(w) = {2na 2 )~ l l 2 e~'i w l a dw. In 
the original coordinate vq for the velocity of the wall, the distribution is 

where := m§o\. 

-v 2 




Figure 3.3: A typical trajectory of the Markov chain for the random scattering process of Figure l3.ll The 
parameters are given in the text. 



The random scattering map given in Definition[5]generates a random dynamical system on 
H!" whose orbits are equivalently described as sample paths of a Markov chain process with 
state space H™. One way to interpret such multi- scattering process is to imagine that a point 
mass m\ undergoes a random flight inside a long channel bounded by two parallel lines (the 
channel walls), these walls having at close range (compared to the distance between the two 
lines) the structure depicted in Figure l3~f1 

Figure [3731 shows a typical sample path of the multi-scattering Markov chain obtained nu- 
merically for the contour function f{z\) = \JR 2 - z\ - \JR 2 - a 2 /4, with a\ = 1, R = 4, and 
masses mo = 80 and mi = 1. The variance is do = 1 and the number of iterations is 10 4 . (In the 
figure we used - U2, so the trajectory is shown in the upper-half plane.) 
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According to Proposition!]] the- stationary probability distribution for u = ( jq, x 2 ) is 

(3.1) dr}{u) = l -=u 2 e~^ u \ l l a2 duidu 2 . 

Expressed in the original velocity variables V\, v 2 of m\, this distribution has the form 

(3.2) dt](v) = ^ m ^ 2 v 2 e-z mM2 d Vl dv 2 = {\[\( m ^f' 2 s 2 e-^" hs2 rfsj (^cosBdej 

where s= \ v\ is the speed of mi and 8 is the angle the velocity of mi makes with the normal 
to the reference plane pointing into the region of interaction. The fact that p is the same in 
both distributions of velocities (for mo and mi) is indicative of (thermal) equilibrium. These 
distributions are illustrated in Figure l3~4l 




Figure 3.4: Factor densities of the stationary distribution dr\ of scattered angle (left) and speed (right) 
as given on the right ot l3.2l above, obtained by numerical simulation of the random billiard 
using a sample run of the Markov chain of length 10 7 . The parameters are mo = mi = 1, 

a\ = 1/2, f{zi) = ^R 2 -z\- ^Jr 2 - a\ /4 (an arc of circle) with a\ = 1 and R = 3. The an- 
alytically derived expressions for the distributions are also shown above (dashed lines) but 
are virtually indistinguishable from those obtained numerically (solid lines). 



4 Differential approximation of the scattering operator 

In this section we denote the scattering operator by P} V indexed by the flatness parameter h, 
and define for all <D in the space of compactly supported bounded functions (H™) 

h 

Other choices of denominator can be more natural or convenient in specific cases, but h in- 
dicates the correct order of magnitude. Our immediate task is to describe a second order dif- 
ferential operator £ to which £ f, converges uniformly when applied to elements of ( H" ! ) . 
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4.1 Definitions, notations, and preliminary remarks 



The notations used below were introduced in Subsection 12.41 and are summarized in Figure 
12.31 In addition, we occasionally use the shorthand ^ v := £/ (£, e). Also, the variable r = (7, c) 
will typically be used to represent the initial position of trajectories, instead of (x, c). Given 
an initial state (r,£ ) with r on the reference plane, let V = V(r, <f) denote the component in 
U™ of the velocity of the return state T(r,f ). Recall that, at the end point, the scattering map 
reflects the velocity back into H" +1 . For trajectories that collide only once with the graph of F, 
it will be convenient to introduce the independent variable x e J" as indicated in Figure 1231 
and use it to express both 7 and V for a given (, instead of writing V(7, £) directly. Note that 

7(x,{)=x+(c-F(x))^, 
whose differential in x is d7 x = I- dF x ® £ v . By a standard determinant formula, 
(4.1) det(dr x ) = l-dF x ($ v ) = l + {n v (x),{ v ). 

Recall that h is the supremum over T" of Igrad^Fl 2 . We wish to study the scattering process 
for small values of h. 

Lemma 2. Let f := v + w e IHP x R k be such that the trajectory with initial state (7, c,() collides 
with the graph ofF only once for all 7 e J n . We regard! as a function of the initial velocity 
and a point x e T", as indicated in Figure UT3\ Let V = V(x, v, w) be the component in Hi 1 of 
the velocity of the billiard trajectoiy as it returns to the reference plane after one iteration of the 
scattering map. Then V = v + 2(\ - 2(2, where 

(i(x, v, w) := (n(x),e) ((n(x), v)e+ {7i(x), w)e- ( v, e)n y (x)) 
(z(x, v, w) := (n(x), v)n y (x) + (7i(x), w)7i y (x) + \ n(x)\ 2 (v,e)e. 

IfF is symmetric, these functions satisfy (\(-x, v, w) = -(i(x, v, w), &{-x, v, w) = (2(x, v, w). 

Proof. This is an entirely straightforward calculation, of which we indicate a few steps. The 
reflection of £ after the single collision with the graph of F at x e J" is naturally given by 
= ^-2(n(x),^)n(x) e H" +1 . This is then reflected by a plane perpendicular to e, resulting in 

77 = £' - 2(<f e)e = f + 2{ n, e) ( (n,J)e- (<f, e)n) - 2 ( {n,J)n + \ n\ 2 {{, e) e) e H^ +1 . 

Now apply the linear projection 77 1-* n Y and use that (£, e) = ( v, e) and (n,() = (H, v) + (n, w) 
to obtain the stated identity relating V and v. For the rest, use (n(-x),e) = (n(x),e) and 
n(— x) = -7i(x). □ 

Lemma 3. Define W(v, h) := - \v\ + \(v, e}\/4\/~h. Then for small enough h (e.g., h < (3/4) 2 ), for 
allxel" and all £ in the set 

D h :={v+ weH'l +1 :\w\< W(v,h)} 

the trajectory with initial vector!; starting at (r(x,£) , c) collides with the graph ofF only once 
and the Jacobian determinant of x >-^T(x,() satisfies det(dr x ) = 1 + (n v (x),£ v ) > 0. 
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Proof. Let £' := (-2(^,n(x))n(x). A sufficient condition for single collision is | > 

\fh. In fact, if there is a second collision elsewhere on the graph of F under this condition, a 
comparison of slopes would indicate the existence of a point where the gradient of F exceeds 
\fh, a contradiction. Using = \i\ and simple algebra, this is equivalent to 

(S',e) 2 >(h/(l + h))\t\ 2 . 

Further elementary manipulations give ,e) = -(v,e) + 2(v,e)\n\ 2 -2(v+ w,~n){n,e). From 
\(n, e) \ < 1, \{v+ w ,~n) \ < (\v\ + \ w\)\n\, and \ n\ < \f~h< 1, we derive 

\($\e)\>\(v,e)\-2Vh{\v\ + \w\)-2h\{v,e)\. 

It follows that 



(l-2h)\{v,e)\-2Vh(\v\ + \w\) >Vhy / \v\ 2 + w 2 

is also a sufficient condition for single collision. Since < \J x 2 + y 2 < \x\ + \y\, yet another 
sufficient condition is 

(4.2) |u/|<(l-2fe)^H-|i;|. 

3\/h 

The inequality l + {n v ,£ v ) > can be rewritten as (n, e)\(v, e)\ > \{H, v+ w)\, which is easily 
seen to be implied by 

\w\< ^ — - -\v\. 

Vh 

But this in turn is implied by inequality 14. 2l for sufficiently small h. For small enough h, we 
may simplify l4.2l by writing the right-hand side as \(v, e)\/4\/~h - \v\. □ 



4.2 The operator approximation argument 

Let $ be a smooth function defined on a subset U c R m . The kth differential d k <t> v of <D at 
v € U is the symmetric fc-linear map on T V U such that d k O v (ei 1 ,...,ei k ) = (Di 1 ---Di k <b)(v), 
where Dj is the directional derivative in the direction of the constant vector field ej. If £ is a 
constant vector field, then 




Let g(s) = <J>(f+ s£ ). In the above notations, the Taylor approximation of g(l) up to degree 2, 
expanded in derivatives of g(s) at 5 = 0, has the form 

(4.3) <P( v + () = 0( v) + dO„(0 + ^d 2 Q> v (i;,S) + 

where \R V (Z)\ = \fo ^-^-(d 3 3>) v+tl: ($,$,$) dt\ < |||d 3 <P|||f| 3 . For the main theorem below, 
where $ will be compactly supported in H™, ||rf 3< l ) || may be taken to be the supremum over v 
of any choice of norm on the 3-linear map at v. 
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We introduce linear maps C : U n+1 U k and A: R n+1 -> U n denned by 
C:= / w* ®wdp(w), A:= / n* (x) ®n(x) d\(x). 

jR k Jj" 

Thus, by definition, Cu = f U k(w, u)wdp( w), and A := f jn A(x) dA(x), where A(x) is the lin- 
ear map A(x)u := (u,n(x))Ti(x). Then A and C are non-negative definite symmetric linear 
maps. 

For convenience of notation, we shall often write below £[•••] = f jn ■■■dA(x). Given a twice 
differentiable function O, let Hess^cD represent the matrix associated via the standard inner 
product (•,•) to the second derivatives quadratic form d 2 Q> v . Let Q : U n+1 -* U q be any self- 
adjoint map and denote A® = QAQ. Observe the identities: 

Tr(,4 Q ) = £• [| Q7T| 2 ] , Tr(A Q Hess v o) = E[d 2 ® v (Qn,Qn)] , (A Q u, u) = e[(QH, u) 2 } 

as well as 

d<S> v A Q u = E[(Qn,u) d<b v Qn\, d 2 $> v {A Q v x , v z ) = E[{Qn, Vi) d 2 <& v (Qn, v 2 )] . 
Similar identities hold for C. In particular, 

Tr(AC) =Tr^ c ' /2 j = J k E[{n,w) 2 ] dp(w). 

The orthogonal projection 771-* rf will be indicated by Q Y :U n+1 -»IR m . If <T> is a function on H"' 
which does not depend on w e U k , then grad y <t> = Q v grad !; <J> and HesSj,<J> = Q Y (Hess„<l>)Q v . 
Similarly, we may define the orthogonal projection Q A to U k . if p. is the Gaussian distribu- 
tion with temperature parameter a 2 (see Subsection l3.ll then C := f R t w* ®wdp(w) = a 2 Q A . 
Observe that A goes to linearly in h. The following assumption is very commonly satisfied: 

Assumption 1. We suppose that the limit A := lim^o A/ h exists. 

Theorem 1. Letp be a probability measure onU k with mean 0, finite second moments given by 
the matrix C and finite third moments. Under Assumption^ define the differential operator 

2 

(£<P)(v) = -4(Agrad„<D, v) + - \{Av, v)+Tr(CA) -Tr(A)(v,e) 2 l (grad„fl>, e) + 

(v, e) 

2( v,e)[( v, e)Tr ( AHess y <D) - 2(AHess y O e, v}] + 2 ( ( Av, v) + Tr (CA) ) (HesSyCp e , e). 
on smooth functions <S>. Then 

lim — = £<D 

h 

uniformly on H'", for each <D e C^°(H'"). When p is the Gaussian distribution on U k with 
temperature parameter a 2 (see Subsection \3.1\ l, then C = a 2 Q A . 

Proof. Recall that the translation invariant probability measure on T" is here denoted by A. 
With the notations of Lemma[3]in mind we define T>%(v) := { w e R fc : \ w \ < W(v, h)}. Then for 

<PeC °°(H™), 

(P h 9)(v)= f f ®(V(r,v,w))dA(r)dp(w) = h+I 2 , 
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where for I\ the integration in w is over T>f,(v), and for I 2 the integration is over D^(y). 
Notice that | 1 < (1 -/i(D^( v))) || <t> || 00 goes to as h approaches 0. 

We now concentrate on I\ . Using Lemma |2] and the form of the Jacobian determinant 
det(dT^) given in Lemma[3] 

<$(v + 2(i -2( 2 ) (1 + S(x, v, w)) dX(x) did(w), 

where S(x, v, w) ■■= (n(x), v + w)/(n(x),e)(v,e) and d = d(x, v, w). To simplify the notation 
we write I\ = J rD> ^ I\(v, w) d^( w), where 

h(v,w) :=E[Q>( v + 2( 1 -2( 2 )(l + 8)], 
and E, defined earlier, indicates average over x. We now use the symmetries: 

(l(-X, V, w) = ~(l(x, V, w), ( 2 (-X, Vy uf) = ( 2 (x, v, w), S(-x, v, w) = -8(x, v, w) 
to write 

0(i; + 2fi-2f 2 )+0(y-2fi-2f 2 ) <P(v + 2Ci -2( 2 )-®(v-2( 1 -2( 2 ) J 

1 o . 

2 2 

Notice that (j are of the order 0(h) in h for each v and w . Each<I>(v+ri) maybe approximated 
by a Taylor polynomial at v up to degree 2 14.31 , 

1 „ 

®(v + T]) = <D(y) + d<J> v Ti+ -d d O v (ri,r]) +R v (t]), 

where |i?„(77)| < | ||<i 3 <t>|| ||??|| 3 - The sum of all terms inside £[•••] has second degree Taylor poly- 
nomial of the form 

P 2 (v,(i,( 2 ) = 0(1;) +2(dQ„(-( 2 + d 2 ® v (( 1 ,(i) + d 2 ® v (( 2 ,( 2 ) -2d 2 ® v (( lt ( 2 )8). 

Keeping only terms in Ii ( v, w) up to first degree in h yields 

h(v, w) =®(v)+2E[d® v (-( 2 +(i8)]+2E[d 2 ® v (( 1 ,(i)]+Ermr( v,w,h), 

where the error term is bounded by a product, |Error| < C®p-$(\v\, \w\)h 3 / 2 ; here C$ is a con- 
stant depending only on the derivatives of <P up to third order and p% is a polynomial in | v\, \ w\ 
of degree at most 3 that does not depend on <J> and h. The linear term in d contributes to 
Ii(v,w) the expression 

2 

-4E [In, v)d<S> v n y ] +- rE[(n, v) 2 + (n, w) 2 - \n\ 2 (v, e) 2 ] d(S> v e+ 

(v,e) 

4 _ _ _ _ 

-E[(n, w)(v,e)d<$ v n y -(n,w)(n, v)dO v e] . 

(v,e) 

Since the measure /i is assumed to have mean (and finite second and third moments), the 
last term above (in which w appears linearly) vanishes after integration over v). There- 
fore, the zeroth and first order terms (in <I>) contribution to h are 

h = ot h \<S>(v)-AE [{n, v)d® v n y ] + 2 E[(n, v) 2 + aJ^(C h ~n,Ji) - \n^(v, e) 2 ] d<t> v e\ + ■■■ 



h(v, w) = E 
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where := p{T) h {v)) goes to 1 and Q, := J, Dj / y x w* ® wdp(w) goes to C as fa approaches 0. 

We now proceed to the second order terms. A similar kind of analysis, where we disregard 
first order terms in w and drop terms in fa of power 3/2 or greater into the error term (this 
involves approximating an overall multiplicative factor (n, e) 2 by 1), yields the second order 
(in <J>) contribution to I\ given by the sum a\ + a 2 , where (separately, so as to fit in one line) 

a x = 2a h {{v, e) 2 E [d 2 ® v (n y , n y )] -2( v, e)E[(n, v)d 2 ® v (n y , e)] } 
a 2 = 2a h E[(n, v) 2 ] d 2 <& v (e, e) + 2E[(C h n,n)] d 2 ® v (e, e). 

Collecting all terms, and using the identities listed for A and C noted prior to the statement of 
the theorem, yields 

h = 0( v) - 4 (grad„<P, Av) + [{Av,v) + a^Tc (C h A) - TrA{ v, e) 2 ] (grad^, e) 

+ 2{v,e)[{ v, e)Tr ( Ao Hess y O) - 2 (HesSyOM e)] 
+2 (( Av, v) + a^Tr (C h A)) (HesSyO e, e) + Error( v, fa) 

where the error term is of order fa 3 / 2 . We have used that the third moment of p is finite to 
ensure that the error term is finite. If <t> e C£°(H™), it follows that as fa -> 0, the quantity 
(I\-Q(v))/h has the same limit as ((Pf,Q>)(v) -Q(v)) / h, which is (£<!>) (v), the convergence 
is uniform, and the limit is (£<£)( v) as claimed. □ 

Recall that H"" 1 " 1 = x U is the decomposition of velocity space into "observable" and 
"hidden" components, with respective projections Q Y and Q A defined earlier. LetA y = Q y AQ y 
and A A = Q A AQ A .We make now an additional but very natural assumption, which holds in all 
the examples discussed in this paper, that A is adapted, according to the following definition. 

Definition 7. The linear map A is adapted if A = A y + A A , in which case a similar decomposi- 
tion holds for A under Assumption^ and we say that A is also adapted. 

For adapted A and for C and a 2 as described at the end of TheoremfT] a 2 = Tr(CA)/Tr(A A ). 
Also recall the stationary measure dr)(v) = p(v) dV(v) described in Proposition [1] whose 
density is g( v) = cv m exp (-j\v\ 2 /a 2 ), where c is a constant of normalization. 

Corollaryl. Let the same assumptions of Theorem\J\hold. Further suppose that k > 1 andthat 
A is adapted. Let e\,..., e m -\,e = e m e R m be an orthonormal basis of eigenvectors ofA y , with 
A y ei = Aei = A ; - e,-, and A m = 0. The partial derivative of a function O on in the direction e ; - 
is denoted <£> i and the coordinate functions are vi := (v,ei). Then, for® eC£°(H r "), 

(Uq>) ( v) = A, v 2 + Tr(CA)) [(— - ^) 9 m (v) + * mm (i/)l + "Z h ( v) 



where Li is defined by 

m-l 



(£,i®)(v) = -2v i ®i(v) + vf n ®ii(v)-2v i v m ® im (v) 



l-^TrCA-))" 1 Z^jvj 

7=1 



V m Q? m {v). 
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This rather cumbersome expression can be greatly simplified by the following coordinate change: 
Xj := Vj for i = 1, ...,m- 1 andx m := \v\ 2 /2a 2 . Leth(x) = 2a 2 x m -x\ x m-i- Then 



dm) (x) = "fit (h(x)O l ) i + M^l e ^ (h(x)e- x ^ 
\2 I i=1 a 



0, 



where <I> is a compactly supported function on { x : 2a 2 x m > x\ h h x } . 

Proof. This is derived from TheoremQ]by straightforward calculations. □ 

As a special case, suppose that n = 0. Then m = 1 and H" J = (-oo,0), in the direction of the 
single vector e. Write A = A > and C = a 2 > 0. Here, only the speed, i/ e (0, oo), is of interest. 
We denote by cp' and <I>" the first and second derivatives with respect to v. Then 

Corollary 2 (Dimension 1). Under the assumptions of Theorem[l\and that n = m = k = 1, then 
for any compactly supported smooth function <I> on (0,oo), 



\ f (7 / 



(4.4) (£0>)(y)=2A(/ 
r/j/s can foe written in Sturm-Liouville form as 

1 1 d / £) d<1 M 

2 Act 2 p \ / 



a 2 v exp 



( _ 2ct 2 )- 



Proof. This is a straightforward consequence of Corollary[T] Note that the coordinates f,- are 
absent for i = 1,..., m- 1 and Tr(CA A ) = Act 2 . □ 

Consider now the case k = 0, or m = n+ 1. This means that only the initial position in T" is 
random, while the initial velocity is fully specified. Then, as the speed | v\ of the billiard trajec- 
tory does not change after collision, we may restrict the state space of the Markov operator P 
to the hemisphere of radius p := \ v\ in H" +1 . This hemisphere is diffeomorphic to the ball Dp 
of radius p, via the linear projection Q : ->■ U n taking e to and fixing the other coordi- 
nate vectors. In this special case, we can restrict attention to functions of the form O = f o Q, 
where 4>( "v) is a smooth function on Dp and D = Q v. For these functions, (grad^O, e) = and 
(HesSy® u\, Uz) = if either ui or «2 or both are multiples of e. Thus the operator £ reduces 
to 

(£ V P)(77) = -4( Q grad-T-, AT7) + 2 (p 2 - |T7| 2 ) Tr( AHess^). 

Corollary 3 (Constant speed). Let the same assumptions of Theorem[l\hold, and that k = 0. 
Without loss of generality, let the particle speed be 1. Let X\ > 0, i = 1, . . . , n and ei be as in 
Corollary^ while i>( is now used as the coordinates on D" whose coordinate vector fields are 
the e,- . In this new system the operator L has the Sturm-Liouville form 

(4.5) (£¥)(!/) = z£ Aiftl-M 2 )^),- 

(=1 
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where the index in indicates partial derivative in V{. In dimension n = I, L is the standard 
Legendre's differential operator on the interval [-1,1] up to a multiplicative constant. 

Proof. This readily follows from the general form of the operator. □ 

When k > 1, define the inner product 

(0,^):= f ^(v)^{v)p(v)dV(v) 

JM m 

on compactly supported smooth functions. When k = 0, we restrict the functions to the unit 
hemisphere equipped with the measure (given by a scalar multiple of) ( v, e) da>( v), where a> 
is the Euclidean volume measure on the hemisphere, and define the inner product accord- 
ingly. (In this latter case, the density of the measure is proportional to the cosine of the angle 
between the vector v and the unit normal to the boundary of H™ .) We say that £ is symmetric 
if(£<5, v P) = (<I>,£ v P). 

Theorem 2. Under the general conditions of Theorem^ assume further that A is adapted and 
positive definite. (Recall that A is in general non-negative definite.) Then L is a second order, 
symmetric, elliptic operator on (H" 1 ) . 

Proof. The claims are obtained by a long but completely straightforward calculation. We only 
check ellipticity for k > 1. (The case k = is even simpler.) Recall that the symbol of the 
second order operator £ is the quadratic form CT£(£) = T, a ij( v )ii^j> where the ajj(v) are 
the coefficients of the second order terms of £ and £ is a vector of dimension m. Starting from 
the expression of £ given in Corollary[T] the symbol can the written in the form 

m— 1 

ffx(f)=2 £ Ai(v m ^-ViU) 2 + 2a 2 Ti(A A )e m . 
i=i 

Since A,- > for i = 1, . . , , m, and both o 2 > and v m > 0, then o £ (£) = only iff = 0. □ 

That L is symmetric and elliptic can be seen more easily by noting that it can be put in 
Sturm-Liouville form relative to the MB -distribution g. To see this, we first introduce the 
following first order differential operators in U m . (The subindex m in <D m and (;-)m indicates 
derivative in the direction e = e m .) For a smooth function cp, 

(D«) ( v) ■= Vz [ A 1 ' 2 ( v m grad„ 0> - * m ( v) v) + Tr ( CA) 1/2 <P m ( v) e] . 

If 3 is a vector field in IR m , 

(D'H)(iO := s/2 [-div(v m A 1 / 2 -) + ( v,A l l 2 S.) m -Tr(CA) 1 ^ 2 <E,e) m ] . 

Then D' is the adjoint of 2) relative to the Lebesgue measure on U m . That is, if either <J> or H is 
compactly supported, then 

f (DQ>)ZdV= f 0>(D'H)dV. 
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We now restrict these operators to the half-space H'" and define D*H := g lr D' (pH) . Clearly, 
D* is the adjoint of D with respect to the density g: 

f (1>*E)(S>gdV= f Z-(T)®)gdV 

JH m JU m 

Proposition 3. Under the assumptions of Theorem^ and that A is adapted, the differential 
operator L has the form 

where® is a smooth, compactly supported function ml™. 

Proof. This amounts to a tedious but entirely straightforward exercise. □ 

5 Diffusion limits of the iterated scattering chains 

One reason for relating the Markov operator P to an elliptic second order differential operator 
is the desire to obtain diffusion approximations of Markov chains associated to our random 
mechanical models. In this section we turn to such approximations. 

5.1 Generalities about diffusion limits 

The results stated here are corollaries of Theorem []] and general facts about diffusion limits 
from [TBI. Chapter 11. 

Let K be an open connected subset of U m . We shall soon specialize to K = H r " after review- 
ing some background information. Let O be the space of continuous functions from [0, oo) to 
K. Define n t :Q.^JC such that n t (a)) =u>(t). Then Q has a natural metric topology making 
it a Polish space, relative to which these position maps are continuous. Let M be the Borel o~- 
algebra on fi, which is also the cr-algebra on Q generated by all the n t . Let Mf be the cr-algebra 
on Q generated by the n s such that Q< s< t. 

Now consider a (time independent) second order elliptic differential operator L with con- 
tinuous coefficients acting on compactly supported smooth functions on K. After [15], a 
probability measure P on (Q,M) is said to be a solution to the martingale problem for £ 
starting from (s, v) e [0, oo) x K if the P-probability of the set of paths w such that u>{ i) = v 
for < t < s is 1 and <p o n t - f (Lq))o n r dx is a P-martingale after time 5 for all compactly 
supported smooth q> on IK. 

Lemma 4. A sufficient condition for the martingale problem to have exactly one solution is the 
existence of(i) a non-negative function <p e C 2 ( 9£ ) such that(p{ u n ) -> oo as u n -* oo (that is, u n 
eventually leaves every compact set as n -> ooj and (ii) a constant A > such that Lip < X(p. 

Proof. The proof is easily extracted from the proof of Theorem 10.2.1, p. 254, of [HI. □ 

Given a family of transition probabilities kernels u i-> LI^ ( u, ■ ) with state space IK parametrized 
by h, define for each v e K the family p'j of probability measures on O characterized by the 
following properties ([15], p. 267): 

1. The set n~^{{v}) has f h v -probability 1; 
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2. The set of polygonal paths oj such that 



w(t) 



(k+\)h-t 



cu(kh) + 



t-kh 



co((k + l)h) 



h 



h 



for all integer k > 0, has P „ -probability 1; 
3. The conditional probability given equals Ylh(ti)(kh),-); that is. 



pS(»(fc+i)fc6r|M fch ) = n fc (» fcfcf r) 



for all fc > and all F in the Borel cr-algebra of "K. 

Conditions 1 and 2 mean that the distribution of (7io,7ii 1 ,n 2 h,- ■ ■ ) is the time-homogeneous 
Markov chain starting from v with transition probabilities u t-» II/;(w,-). Notice that we have 
used before the notation v u for n^(w,-). Let Pf, be the corresponding operator on compactly 
supported smooth functions and let Ah ■= Ph ~ I- Condition 3 above is equivalent to: 



for every compactly supported smooth function cp on Oi. 

The key fact we need from the general theory of diffusion processes can now be stated. 

Theorem 3. Assume that (i) the elliptic second order differential operator L (with continuous 
coefficients) is such that for each u e "K there is a unique solution P u to the martingale prob- 
lem for £ startingat u; and (ii) h~ x A] % <& converges to £<t> uniformly on compact sets for every 
smooth compactly supported <J> on "K. Then lim; 2 ^ = and convergence is uniform in u 
over compact subsets ofK. 

Proof. A proof is easily adapted from that of Theorem 11.2.3 of fl5l. □ 



We now set 'K = IHl!! ! . It was shown above that for the convergence of the Markov chain to a 
diffusion process it is sufficient to have: (i) convergence of h A},® to £<t> for every compactly 
supported smooth <J> as in Theorem|4]and (ii) a function q> as in Lemma|U The convergence 
required in (i) is implied by Theorem[T] We now show the existence of a (p. 

Lemma 5. Let L be the differential operator of Theorem^ Suppose that A is adapted. Then 
there is a smooth function <p : H™ -* (0, oo) and a positive constant A such that Lq) < Aq> and 
q>( v) -* oo as | v\ oo orv approaches the boundary of¥\™. 

Proof. As A is adapted, we may assume that L is as in Corollary [T](fc > 1) or as in Corollary|3] 
{k = 0). The case k = is much simpler: take (p(v) = c + \ v\ 2 for a big enough constant c. So we 
assume £ is as in Corollary[T] Let u : (-oo,0) -*■ (0, 1] be a smooth function such that u(s) = 1 
for \s\ > 1 and u(s) = -s for \s\ < 1/2. Now define 



fc-i 

V ° Kkh -T t ( A h ( P)° 7l jhisa (M kh , Pj5)-martingale 

7=0 



5.2 Back to the random scattering operators 



<p(v) = c + 8Tr( CA) + 2Tr( A) + 1 v | 2 - In u{ v m ) 
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where c is a positive constant still to be chosen. It is clear that q>( v) goes to infinity as claimed. 
A straightforward computation shows 




8Tr(CA)+2Tr(A) if|i/ m |<l/2 
8Tr(CA) ifK|>l 



If 1/2 < \v m \ < 1, the coefficients of Lq> are seen to depend quadratically on V\,..., v m -\ and 
are bounded in v m ,\\ v m . This shows that in this range of v m there are constants c, X greater 
than 1 such that Lcp < c + A| v\ 2 . On the other ranges, L<p < <p, in which c = 0. □ 

Thus we conclude: 

Theorem 4. The martingale problem for the random billiard differential operator L under 
the assumption that A is adapted has a unique solution = P> u for each u e H™. Furthermore, 
lim^o = Put where P„ solves the martingale problem for the Markov chain with transition 
probabilities operator Pf, . Convergence is uniform in u over compact subsets ofK. 

Proof. This follows from Theorem[3]and Lemmas[4]and[5] □ 

It is useful to express the diffusion process with infinitesimal generator L as an stochastic 
differential equation. 

Proposition 4 (Ito SDE). We consider separately the cases k=Q and k > 0. The operator A is 
assumed positive definite on K m_1 . 

1. Under the conditions of Corollary\^(k = 0), the ltd differential equation associated to the 
infinitesimal generator L has the form 

d V t = -4 AV t d t + [2 ( 1 - 1 V t | 2 ) A] 1/2 dB t , 

where B r is n -dimensional Br ownian motion restricted to the disc D". The Lebesgue mea- 
sure on the disc is stationary for this process. 

2. Under the conditions of CorollaryUl (k > 0), the same ltd differential equation has the 
form 

dV t = Z(V t )dt + b(V t )dB t , 
where B t is m-dimensional Brownian motion restricted to H™ , Z( v) is the vector field 

Z(v) := -2Av + (—- ^] ({Av, v)+Tr(CA))e m 
\v m a L I 

and b(v) is the linear map 

b( v) u ■= v m A 1/2 u - ( A 1/2 v, u) e m + Tr ( CA) 1/2 u m e m . 
The stationary distribution is Maxwell-Boltzmann, as described in Proposition^ 
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Proof. This is an easy exercise. The general relation between the infinitesimal generator of 
the diffusion and the Ito equation can be found, for example, in [14] . When k = 0, note that 
the drift term always points into the disc as (Ay, v) > for all non-zero v. As L is a symmet- 
ric operator (see Theorem[2), J (£<£)( y) d[i(v) = for all compactly supported smooth <p, 
where \i is the stationary measure for P given in PropositionQ] The claim about the stationary 
distributions is a consequence of this property. □ 

Proposition 5. Let L be as in Theorem^ Then the following two conditions are equivalent: 

1. L\v\ 2 = 

2. Tr(CA) = and Tr(A) = Tr(A Y ). 

If these conditions hold, the diffusion associated to L restricts to hemispheres of arbitrary radius 
(i.e., the level surfaces of\ v\ z ), and is equivalent to a Legendre diffusion. 

Proof. This follows from the observation that L\v\ 2 /4 = 2Tr( CA) + v m (Tr ( A Y ) - Tr( A) ) . □ 

The significance of this remark is the following. In the examples, Tr(A) -Tr(A v ) consists of 
mass ratios whose denominators are the masses associated to the velocity covariance matrix 
C. These constitute the "wall subsystem," whose kinetic variables are "hidden." Therefore, the 
two conditions amount to the assumption that the masses of the wall subsystem are infinite 
and have zero velocity. Thus we have an elastic random scattering system. 

5.3 Example: wall with particle structure 

Consider the idealized physical model depicted in Figure I5T1 It consists of k point masses 
mi,..., mjc that can slide without friction on the interval [0, 1] independently of each other, 
and a point mass m that can similarly move on the interval [0, oo). When reaching the end- 
points of [0, /], masses m,- bounce off elastically, while m collides elastically with the m, but 
moves freely past z = I. One may think of the m, are being tethered to the left wall by imagi- 
nary (inelastic, massless and fully flexible) strings of length Z; when a string is fully extended, 
the corresponding mass bounces back as if due to a wall at z = I. So the m; are restricted to 
[0, /], but m is free to move into this interval and may collide with the m,-. 







I 



Figure 5. 1 : In this model, masses m; constitute the wall system, while m is the free mass. The system is 
essentially one-dimensional. 
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The positions of the m; are Zj e [0, Z] and z e [0, oo). Let M = m + mi H H ra^. In the new 

coordinates x,- = \pm\ fM z,-, X£ +1 = \J~m\Mz, the kinetic energy form becomes 

= (M/2)(*? + - + 4 Hl ). 

We may equivalently assume that {x\, . . . , x^) defines a point on the torus J k by taking the 
range of x\ to be [-at/2, at/2], where a, = 2\J mi/Ml, and identifying the end points a,/2 and 
-ail 2. Mass m is then constrained to move on the interval defined by 

x fc+ i > F(xi, . ..,x k ) ■■= maxjv' m/mi \xi\, . . . ,y/ m/ m k \x k \} . 

Thus the configuration manifold is M = { (x, x k+i ) e T x U : x^+i >F(x)\, and collision is rep- 
resented (due to energy and momentum conservation and time-reversibility), by specular re- 
flection at the boundary of M as depicted in Figure lfT2l 

This deterministic billiard system can be turned into a random scattering system in several 
ways. We illustrate two natural possibilities, which we call the heat bath model and the ran- 
dom elastic collision model. The assumptions for the heat bath model are as follows: at the 
moment m crosses z = I into [0, 1], the initial position of each m; is a random variable uni- 
formly distributed over [0, 1], and the velocity of m; is normally distributed with mean and 
variance a\. We assume that these aj are such that m^cr? = mja 2 j for all i and j. In physical 
terms, we are imposing a condition of equipartition of energy among the wall-bound masses. 
In the new coordinates x\ the velocities are normal random variables with mean and equal 
variance a 2 = (mt/ M)a 2 . 




This is essentially the case described in CorollaryFJ] The dimensions for the heat bath model 
are: n = k = m- 1. Let e\,...,ei c ,e = e^+i be the orthonormal basis of coordinate vector fields 
corresponding to the x ; - . Then the projection to the hyperplane e 1 of the normal vector field 
n(x) to the graph of F is 



n(x) = ± A / e\ for x such that > max|x, |. 

V m+mi j 
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Therefore, A is the diagonal matrix 



A=V Y I 



m + m: 



where Vq is the volume of the sector |x, | > maxy \xj |, normalized so that the total volume of the 
torus is 1. The matrix C is the covariance matrix of the velocity x, and is by assumption the 
scalar matrix C = a 2 e* ® e,. 
Therefore, 

Tr (CA)/Tr(A) = a 2 , Tr( A) = V" , h = max — , . . . , [. 

f^m + rrii j [rm m k ) 

Let us say for concreteness that all the wall-bound masses are equal to mo, so h = m/ mo and A 
becomes the identity matrix times Vq, whose trace is kVo. From Corollary[T]we obtain, for the 
heat bath model with small ratio m/mo and k equal bound masses, the following differential 
operator. Let v indicate the velocity of the free mass m (in the new coordinate system, so 
v= Xfc+i) and let <p be any compactly supported smooth function on the interval (0,oo). Then 

(5.1) Lt> = 2kV a 2 

The corresponding Ito diffusion has the form 

dV t = 2kV a 2 (---^-) dt + \/2kV a 2 dB t . 
\v a 1 ) 

Figure [573] shows a sample path for this SDE obtained by Euler approximation. (See Q2] ■) 

We now consider the random elastic collision model. The assumptions for this model are 
as follows: the velocities of all the masses m,-, m constitute the observable variables, and the 
positions in [0, 1] of the wall-bound masses at the moment m crosses into [0, Z] are uniformly 
distributed random variables. This is the case to which Corollary|3]applies, where n = k. (The 
integer k of Theorem[T]is 0.) Again, for concreteness, suppose that all the wall-bound masses 
are equal to mo. The eigenvalues of A = Vol are all A,- = Vq. We may assume without loss 
of generality that the constant speed of the billiard particle (in the billiard representation of 
Figure [5721 is 1 and let v denote the projection of the billiard particle's velocity to the unit disc 
D\ in dimension k. Then we obtain from Corollary|3] 

k 

(5.2) (£0)(i/) = 2VbX((l-M 2 )<I>.-)j 

(=1 

where <D is any compactly supported smooth function on D\ . 

The differential operator of 15.21 as well as 14.51 in Corollary [3] generalize in a natural way 
the standard Legendre operator defined on functions of the interval [-1,1]. We refer to the 
diffusion process with this type of infinitesimal generator a (generalized) Legendre diffusion. 
A sample path is shown in Figure [T73l 



30 



It is interesting to notice that the heat bath and random elastic collision models lead to very 
standard Sturm-Liouville differential operators. For the heat bath, Equation l5.1l is, up to con- 
stant, Laguerre's differential operator. Essentially the same model of heat bath/thermostat 
described here is used in |4j to build a minimalist mathematical model of a heat engine, de- 
scribed as a random system of billiard type. 



Figure 5.3: A sample path of the Ito equation dVt = (1/ v - v)dt + dBt- We have used Euler approxima- 
tion with time length 10, initial point v = 10, and number of steps = 10000. The mean value 
relative to the stationary distribution is \Jn\2, corresponding to the horizontal dashed line. 



5.4 Example: collisions of point mass and moving surface 

Here we give the differential equation approximation of P - I for the example of Subsection 
13.21 The main interest in this example is that it is the simplest that combines the features of 
the two cases considered above in Subsection l5.3l 

It is first necessary to describe the operator A (see Subsection l4.2l . The notations are as in 
that subsection. The torus T 2 has fundamental domain (centered at (0,0)) 

<20 / TUn 

\x \<t/2, jci < 1/2, t=— W— - 
a\ Y m\ 

The billiard boundary surface is the graph of xz = F(xo,Xi), where 

F(xo,xi) = /— |xo| + a^ l f(a 1 xi). 
V m o 

Let eo, ei,e:= ez be the standard coordinate vector fields for the coordinate system (xo,xi,X2). 
It is easily checked that A is 

A=i^ + 0{h 4 ) S je* Q ®eQ + ^£[f'{a l s)fds+0{h i )^e* l ®e l , 
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where the error term satisfies < 0(h 2 ) < h 2 , while the norm of A satisfies ||j4|| < h. As an 
example, take 

f(zi) = ^R 2 -Z 2 -^R 2 -a 2 /4. 

The graph of / is an arc of circle of radius R intersecting the Zi-axis at the points (±ai/2,0). 
Let the scale-free curvaturebe k := a\/R < 1. Then 

a:= f o \f'(a 1 s)] 2 ds = K-Hn^-l= 1 ^ + 0(K 3 ). 

Thus for small values of h (disregarding terms in k or order greater than 2, and in mil ma of 
order greater than 1) we have 

A = — e ® e + — e x ® e\ and h = — + — . 
mg 12 4 mo 

The operator C takes the form 

/oo 
w 2 dfi( w)eQ ® eo = o 2 eo ® en. 
oo 

We observe that Tr( CA) = j^o~ 2 and Tr( A) = j^ + a. For the special case of an arc of circle, a = 
?c 2 /12, where k is the scale free curvature. This yields the approximation, written informally 
as 

(P<S>)(v)-<S>(v) mi 2 k 2 

(5.3) ^ ^ ^ * - J.^ $ + — 

2 mo 12 

where 

£,em P <P = ( — - ^ ) + 22 



\V2 a 2 ) 



y2 _ y2 

^curv® = -2V\<&1 + — -$2 -2f 1 l'2 , l ) 12 + ff^ll + 

^2 

The mass-ratio and curvature parameters may a priori go to independently with h (under 
Assumption!]] and the particular way in which each goes to matters for the limit. Expression 
!5.3l shows that taking h for the denominator in the quotient (Pf,<S>-<5?) / h used in the definition 
of £ is essentially an arbitrary choice. One could have taken instead Tr(A), for example. If we 

further ask in this example that the scale-free curvature and the mass ratio be coupled by a 

2 

linear relation such as ^ = a^-, for a fixed but arbitrary constant a > 0, and keep the original 
choice of denominator h, then 

a 1 

A= lim A/h= en ®eo + — ; r^i*® e\. 

h-*o 1 l + a 3(1 + a) 



This gives the family of operators (depending on a) 

w \ (P h O)(v)-0(v) 2a 2 a(/l v 2 \ 

(£0) v) = lim V A / ^ = § O2 + $22 

fc-o h l + a l\ V2 a 2 1 



v 2 



2vi®i + — -<P 2 -2vi v 2 ®i2 + v\<S> n + v 2 ®22 



3(1 + a) { v 2 
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In the concrete example of Figure [L2l we took a = 1, a 2 = 1/3 (and multiply the the operator 
by an overall factor 3 to make it look simpler). 

The expression !5.3l points to a separation between, on the one hand, the term responsible 
for the change in speed, which contains the variance (temperature) a 2 and the mass ratio, 
and on the other, a purely geometric term that involves the scale free curvature k. If we let 
a 2 be and the wall mass oo, and consider k to be small, the MB-Laplacian reduces to £ curv . 
It is interesting to note that £ cu „|i/| 2 = 0, so the diffusion associated to this second order op- 
erator restricts to hemispheres of arbitrary radius, and we have a Legendre diffusion. (See 
Proposition[5]) 

Acknowledgment: Hong-Kun Zhang was partially funded by NSF Grant DMS-090144 and 
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